Four-wave mixing of linear waves and solitons in fibres with higher order dispersion 
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We derive phase-matching conditions for four-wave mixing between solitons and linear waves 
in optical fibres with arbitrary dispersion and demonstrate resonant excitation of new spectral 
components via this process. 



The use of any optical system, including fibres, for fre- 
quency conversion relies on the ability to satisfy so-called 
phase matching conditions, which critically depend on 
the dispersive properties of the system, and on strong 
enough nonlinear properties, to allow reduction in the 
threshold pump power. The recent surge of interest in 
theoretical and experimental studies of optical paramet- 
ric processes in photonic crystal fibres (PCFs) [1-7] is 
related to their high nonlinearities, achieved by reduc- 
tion in the core size [8], and to the possibility of disper- 
sion control by suitable design of the photonic crystal 
cladding [2]. 

Here we study parametric generation of new frequen- 
cies resulting from four-wave mixing (FWM) of solitons 
and continuous waves (cw's) in optical fibres under condi- 
tions where the effects of higher order dispersion are im- 
portant, i.e., when pulses are short and/or the frequency 
dependence of the group velocity dispersion (GVD) is 
steep. Interest to this problem arises from the fact 
that the Fourier components of a soliton are dispersion- 
less, while freely propagating cw's are strongly disper- 
sive. Therefore the phase-matching conditions are ex- 
pected to be satisfied at frequencies different from those 
generated by the mixing of cw's. Addressing this issue 
is timely because of the availability of strongly nonlin- 
ear small-core PCFs [7,8], which decrease the thresh- 
old for observing of parametric processes by one to two 
orders of magnitude compared to conventional fibres: 

Iconven ~ lO^W" 1 ™.- 1 , 7 pc/ ~ 10" 2 to lO"^" 1 !^ 1 , 

where 7 is the nonlinear fibre parameter [9]. Strongly 
nonlinear PCFs have already been used to demonstrate 
the coupling of solitons and cw radiation in: supercon- 
tinuum generation [1,6], strong red and blue resonant, or 
Cherenkov, radiation from solitons [6,7] and cancellation 
of the soliton self- frequency shift by the spectral recoil [7] . 
Note, that it is natural to expect that FWM of solitons 
and cw radiation is one of the many nonlinear processes 
contributing to the shape of the supercontinuum spectra. 

The problem of mixing of solitons and cw using the ide- 
alized nonlinear Schrodinger equation (NLS), i.e., with 
higher order dispersion disregarded, has been analyzed 
in a number of papers in the past [10-15]. In these cases 
several exact analytical solutions for a soliton sitting on 
the cw background were found [14,15] and different per- 
turbation techniques suggested [10-13]. However, none 
of these studies have addressed the issue of generation of 



new frequencies by FWM of solitons and cw light - the 
central focus of this paper. 

We assume that the dynamics of the dimensionless am- 
plitude A(t, z) of the fundamental fiber mode is governed 
by the generalized NLS equation [2] 

f+ca 

d z A = iD(id t )A + iA R(t')\A(t~t',z)\ 2 dt'. (1) 
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The dispersion operator in Eq. (1) is given by 
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where r is the pulse duration and luq is the reference fre- 
quency. To avoid any ambiguity in the analytical expres- 
sions we adopt the convention of using round brackets (..) 
to indicate the arguments of functions or operators and 
[..], {..} for all other purposes. R(t) is the response func- 
tion of the material, which includes instantaneous Kerr 
and delayed Raman nonlinearities: 



R(t) = [1 - 6}A(t) + BaQ{t)e' t/T2 sin(Vn) 



(3) 



Here A(f) and 8(t) are, respectively, delta and Heavi- 
sidc functions, a = [t\/t% + T2/Yi]/t2, — 0.18, t\ — 
12.2/s/r, r 2 = 32/s/t [9]. t is the time in the reference 
frame moving with the group velocity vq — v(u>o) and 
measured in the units of t: t = [T — z/vq]/t, where T 
is the physical time, z — Z/L gvc i, where Z is the dis- 
tance along the fiber and L gv d — t 2 /|/32(wo)| is the GVD 
length. The field amplitude A is measured in units of 
N/ J 7 Lg V d, where N 2 is the ratio of the peak power of 
the pump pulse to the peak power of a fundamental soli- 
ton with duration r. 

The dispersive properties of the soliton and linear cw 
are crucial for the following, so we now discuss them in 
some detail. Looking for a linear wave solution of (1) in 
the form A - e iD*z-iS„t we find Dx = rj(S x ). In what 
follows the subscript x can take any convenient notation. 
E.g., S s and S cw correspond below to the frequency shifts 
of the soliton and cw-pump, respectively. Thus the physi- 
cal wave number of a linear wave with frequency ioo+5 x / r 
is given by k x — k{uJo)+D x /L gv d- Plotting k x — fc(ojo) vs 
S x one simply recovers the dispersion profile of the fibre. 
The single soliton solution 

A = F ^) e -iSst+i[D s+q ]z ^ F (q = ^/2q- sech ^/ r ) (4) 
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satisfies Eq. (1) if all derivatives of the function F(£) 
higher than second are disregarded, 8 — 0, and D" < 0, 
i.e., the GVD at the soliton frequency is anomalous. Here 
£ = t — D' s z, r = y/—D"/ (2q), 1 denotes the derivative 
with respect to 8 and q > is the additional shift of 
the soliton wave number. Representing F(£) through 
the inverse transform of its Fourier image F(8), i.e. 
F = J dSFe l ^ Ss ~ s \ one finds that the wave number of a 
Fourier component of the soliton with frequency u)q+8 x /t 
is given by k s/x = k(u> ) + {F> s + q - [S s - 8 X ]D' S }/ L gvd . 
Thus plotting k s / x — k(u>o) vs 8 X will give the tangent line 
to the curve k x — k(u>o) taken at the point 6 X — S s and 
shifted up by q. The linear dependence of k s / x on S x is 
a reflection of the dispersion-free nature of solitons. 
We seek solutions of Eq. (1) in the form 



A = {F(Z)+g(z,Q}e 



iz[D s +q]-iS B t 



£ = t- zD' s . (5) 



Assuming that g is a quasi-linear wave we derive 

ip = d z g - W(id s )g - i2F 2 g - iF 2 g\ (6) 

where D{id{) = —q — D s — iD' s d^ — D(id^ + S s ) and 
p = [D{id£ + S s ) -D„- iD' s di. + l/2D'^d 2 ]F. Now we 
split g into two parts 

g = we i4 > + 1>, <j> = zD(8 cw - 8 S ) + ^[8 S - 8 CW ] (7) 

The ui-term in Eq. (7) is the weak cw pump, which 
obeys Eq. (6) with p = F = 0. The tp-term is the 
field generated through the mixing of the soliton and cw 
pump and by the soliton itself. Seeking ip in the form 

V> = Ipp + ^ +e iz£> ^- S ^ + ^*_ e -izD(5 cw -5 s ) we find that 

ip p obeys 
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and tp± are governed by the system of coupled equations 
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iF 2 we i ^- 5 ' 
where 

t = 

W = D(S CW - S s ) 



1>- 



(9) 



D(id 6 ) - 2F 2 -F 2 

F 2 D{-idz) + 2F 2 
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Eq. (8) does not depend on the cw-pump and it is known 
to describe the emission of non-localized dispersive waves 
from solitons so-called resonant or Chcrcnkov radiation 
in fibres with zero GVD points [6,7,16-18]. This radi- 
ation exists because the continuous part of the spectrum 
of the operator C has a zero eigenvalue, which ensures 
resonance with the forcing term p. The operator C + W 



is expected to have several different resonances, a situ- 
ation which to our knowledge has not previously been 
considered. These resonances are driven by FWM be- 
tween the soliton and the cw, see the left-hand side in 
Eq. (9). To find continuous spectra of C and C + VV we 
neglect the i 7 ' 2 -terms and look for cigenfunctions in the 
form V P ~ e iKZ+t ^~ 5 ^ and ip± ~ e ±i^±i[«.-«]€ Assum- 
ing that the wave numbers of the continuum modes are 
matched with the wave numbers of the driving terms in 
Eqs. (8,9), i.e. k = 0, we find: 

q + D s -[6 s -6}D' s = D(S), (10) 
±D CW T {q + D s - [S s - S CW ]D' S } + 

{ q + D s -[6 s -6]D' s } = D(6). (11) 

Eqs. (10,11) are the equations for 8. The roots of these 
equations, which can be easily found graphically, yield 
the frequencies of the emitted radiation. The right-hand 
side D(5) is simply the dispersion of linear waves, which is 
a nonlinear function of 8. The left-hand sides are straight 
lines. Fitting the experimentally measured dispersion 
characteristics of the fibres usually results in a very high 
order polynomial for D(8) [2,7]. Here, however we use 
the simplest illustrative example, when only the third 
order dispersion is included, i.e. D(S) = —S 2 /2 + e<5 3 , 
e = 8^02 (wo)/(6t|/?2(wo) I)- This model accurately de- 
scribes regions where the frequency dependence of /^(w) 
is quasi-linear - typical for telecom and PCF fibres, see, 
e.g., [7]. Figs. 1(a) and (b) show plots of the right- 
and left-hand sides of Eqs. (10,11). Eq. (10) gives the 
Cherenkov resonance [6,7,16-18], which does not depend 
on the cw-pump. The roots resulting from Eqs. (11), 
however, are determined by the cw-pump. Depending on 
the value of 8 CW Eqs. (11) produce either two or four new 
roots, see Figs. l(a),(b). Assuming that 8 = S cw , we find 
that Eq. (11) with + is satisfied, which means that one of 
the new roots always coincides with 8 CW . Clearly the po- 
sition of these new resonances strongly depends on, and 
can be controlled by, the frequency of the cw-pump and 
by engineering the linear dispersion profile of the fibres. 

The four-wave mixing nature of the new resonances 
becomes clear if we make use of the expressions for the 
wave numbers of the dispersive waves and the Fourier 
components of the soliton (see the discussion around Eq. 
(4)). Indeed, Eqs. (11) can be rewritten in the form 



(12) 



Here k ra d, k cw are wavenumbers of the resonant radiation 
and the cw-pump, and k s / rad , k s / cw are the wavenumbers 
of the Fourier components of the soliton at the frequen- 
cies of the resonant radiation and the cw-pump. 

To confirm our analytical findings we have carried out 
a series of numerical experiments with parameters close 
to the ones in [7]. Comparing Figs. 1 (a),(b) with 1 
(c),(d) one sees excellent agreement between the analyt- 
ical predictions and the frequencies emerging from the 
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modelling of (1). The efficiency of excitation of the new 
frequencies strongly depends on the choice of S s and 5 CW 
relative to each other and to the zero GVD point. This 
explains why not all the resonances are observed simulta- 
neously. Changing 5 cw ,s we have been able to observe all 
of the newly predicted resonances. Theoretical analysis 
of this problem is possible within the framework of Eq. 
(7) and we leave it for a future study. The cw powers 
required to observe the new FWM resonances are of the 
order or less than 1W, see the Fig. 1 caption. 



wavelength x 10" 6 [mj wavelength x 10 6 [m] 

1.480 1.355 1.249 1.159 1.480 1.355 1.249 1.159 




FIG. 1. Resonance frequencies (a,b) calculated from 
Eqs. (10,11) and supporting numerical modelling (c,d). 
Straight /curved lines in (a,b) are the left/right-hand sides 
of Eqs. (10,11). Choosing lj = 2tt x 240THz and the 
fibre parameters from [7] we find /?2 = — 47ps 2 /km and 
d^/3 2 = -0.7ps 3 /km. S s = 0, w = 0.1, S cw = -15 for (a,c) 
and 8s — —4, w — 0.3, 8 cm — —40 for (b,d). Other parameters 
are r = 170fs, q = 4, N = 1.5. For 7 = 0.05(Wm) _1 this cor- 
responds to the pulse peak power 225W, the cw-power 0.12W 
for (c) and 1.2W for (d). The full lines in (c,d) corresponds 
to 8 = 0, i.e. Raman is off, and the dashed lines to 9 = 0.18, 
i.e. Raman is on. The propagation distance is 2.8m (c) and 
16m (d). The top axes are marked in wavelength units. 

We have analyzed FWM between solitons and cw- 
pump in fibres with higher order dispersion and pre- 
dicted the generation of new frequencies, which can be 
controlled by tuning the cw-pump. In addition to their 
fundamental significance, our results have important im- 



plications for the generation of laser light and broad-band 
supercontinua in new frequency ranges. 

This work is partially supported by the EPSRC grant 
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